Black hole formation and slow-roll inflation 
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Black hole formation may occur if the spectrum of the curvature perturbation £ increases strongly 
as the scale decreases. As no such increase is observed on cosmological scales, black hole formation 
requires strongly positive running n of the spectral index n, though the running might only kick 
in below the 'cosmological scales' probed by the CMB anisotropy and galaxy surveys. A concrete 
and well-motivated way of producing this running is through the running mass model of slow roll 
inflation. We obtain a new observational bound n' < 0.026 on the running provided by this model, 
improving an earlier result by a factor two. We also discuss black hole production in more general 
scenarios. We show that the usual conditions e < 1 and \rj\ <C 1 are enough to derive the spectrum 
T-'f(fe), the introduction of higher order parameters £ 2 etc. being optional. 



PACS numbers: 26.35.+C, 98.80.Cq, 98.80.Ft 
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II. FORMING BLACK HOLES 
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The primordial curvature perturbation £ is only of or- 
der 10~ 4 on cosmological scales, but it might be of order 
1 on smaller scales. Primordial black holes will then form 
as those scales enter the horizon, with possibly observ- 
able consequences. The purpose of this paper is to see 
to what extent the value of order 1 is reasonable, taking 
into account observational constraints and current think- 
ing about the origin of (. 

In Section II we see what is required for black hole 
formation, in terms of the spectral index n(\nk) = 
dlnV^/dlnk, which specifies the scalc-dcpcndcncc of the 
spectrum of £. Averaged over the whole range of scales 
we need strongly increasing n corresponding to running 
n' ~ 10~ 2 . Up to this point we assumed nothing about 
the origin of £. In Section III we introduce the assump- 
tion that it originates from the inflaton perturbation dur- 
ing slow-roll inflation (the standard paradigm). Within 
this paradigm, the only extant model giving the required 
running n' is the running mass model, which typically 
makes n' roughly constant hence requiring n' ~ 10~ 2 on 
cosmological scales. 

In Section IV we ask whether such a large value of 
the running is still permitted by current data, thereby 
updating an earlier work. We find that it is. 

The question then arises, whether black hole formation 
can still be achieved if n' is negligible on cosmological 
scales, as might be required by future data. In Sections 
V to VIII we show that black hole formation can indeed 
be achieved within the standard paradigm. Along the 
way, we are led to consider the standard paradigm in 
more detail than before. 

In Section IX we depart from the standard paradigm, 
by allowing a curvaton-type mechanism to contribute to 
the curvature perturbation. We show that black hole 
formation can occur if there is a switch from the standard 
paradigm to a curvaton-type paradigm as we go up in 
scale, or vice-versa. We conclude in Section X. 



A. Viable black hole formation 

The curvature perturbation £ is time-independent dur- 
ing any era when there is a unique relation P(p) between 
pressure and energy density From the success of 
the BBN calculation, we know that such is the case to 
high accuracy a few Hubble times before cosmological 
scales start to come inside the horizon. On cosmological 
scales, the spectrum V((k) is then observed to be about 
(5 x 10" 5 ) 2 H* 1 

When smaller scales start to come inside the hori- 
zon V(;(k) could be bigger. To discuss that case, recall 
that the typical value of £(x) in the observable Universe 
smoothed on the scale A: -1 , is of order V^' 2 {k). If 
Vc^' 2 [k) is bigger than 10~ 2 or so, then black holes will 



form with an abundance that can be ruled out [j| by 
a variety of observations. A somewhat smaller value, say 
V c 1/2 ~ 1(T 3 , would give an abundance whose effect may 
be observable in the future. We want to see how such a 
value may be generated over some range of k. 
The spectral index n is defined by 



n(k) - 1 



d\nV Q {k) dlnV c (N) 



dhik 



dN 



(1) 



In the final expression we assume almost-exponential in- 
flation, with N(k) the number of e-folds of inflation re- 
maining after the epoch of horizon exit k — aH for the 
scale k. We will freely use N as an alternative variable 
to \nk. 



* x The precise number refers to the pivot scale defined in Section 
III. As usual, k is the coordinate wavenumber so that k/a is 
the physical wavenumber, with a the scale factor of the Universe 
normalized to 1 at present. The Hubble parameter is H = a/a 
and horizon entry is defined as k = aH. 
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We take No = 50 unless otherwise stated, where the 
subscript denotes the epoch of horizon exit for the 
present Hubble scale k = Hq. This is the largest observ- 
able scale, and smaller scales leave the horizon AN > 
e-folds later. For scales probed by WMAP and galaxy 
surveys, AN < 7 0. 

Until Section IXB we will assume that the black holes 
form on the scale leaving the horizon at the end of infla- 
tion, corresponding to N — 0. We will take the criterion 
for significant black hole formation to be 7-^(0) = 10~ 3 . 
Then we need hi[V 6 (0) /V C (N )] ^ ln(10- 3 /10" 9 ) ~ 14. 
With constant n this requires n — 1 ~ 14/ No ~ 0.3. This 
was compatible with observation for many years but is 
now excluded. 

Taking instead n' = dn/ dink to be constant black hole 
formation requires 



14 = In 



n(o) 



V C (N ) 



l)+ l -Ntn>. 



(2) 



Since observation requires no — 1 ~ —0.05 Q, the first 
term is negligible and we need n' ~ 28/Nq ~ 0.01. As 
we will see, this is compatible with observation. At the 
end of inflation, n onc i — 1 ~ n' x No ~ 0.5. This might be 
regarded as more or less compatible with the requirement 
\n — 1\ <C 1, that must hold if £ at the end of inflation 
is generated by the perturbation in a single light field 
0,13, HJ. We are going to assume such a scenario and 
therefore rule out |n en d — 1| > 1. 

Finally, suppose that n' increases monotonically as we 
go down in scale. A significant increase would require 
I "•end — 1| <L 1 which we rule out. A significant decrease 
would require n' Q significantly bigger than 0.01, which as 
we will see would be in conflict with observation. 



B. The case P c > 1 

We end this section by discussing briefly the case 
Vc,{k) ^S> 1. If £ is generated from the inflaton perturba- 
tion during slow-roll inflation (the standard paradigm), 
this is ruled out. The reason is that the regime Vq si 1 
then corresponds to eternal inflation [9J , whose duration 
is indefinitely long. Then the slow roll model has noth- 
ing to do with the observed perturbations, which instead 
have to be generated when the eternal inflation is over. 

However, such a perturbation could also be generated 
by the perturbation of a curvaton-like field [1, E(|, as 
one can readily understand from the non-perturbative 
formula 

[E IS Hi C = SN which makes sense no matter 
how big is C. In that local observer would notice 

nothing amiss before horizon entry, and it is not clear 
what will happen at horizon entry. 



such models that have been proposed, the only one giving 
the large positive running required for black hole forma- 
tion is the running mass model [H El El El H3] #2 • 
This model invokes softly broken global supersymmetry 
during inflation, with a potential 



V 



V ± im 2 



(3) 



and a running mass m 2 (</>) whose form is determined by 
Renormalization Group Equations (RGE's). 

Over the limited range of cosmological scales, n(N) 
typically has the two-parameter form 



n(N) - 1 



se c(N -N) _ ^ 



(4) 



With 



c-10" 1 to 10~ 2 (5) 

this gives 

no = 2(s — c) — 1, n' Q = 2sc, (6) 
whose inverse is 

We see that significant negative running is forbidden;^ 3 



n' > -(n - 1)' 



-3 x 10" 



(8) 



Higher derivatives 
pressed; 



t (m) = d m n/d(-N) rn are sup- 



n 



(m+l) 



cn 



(m) 



m > 1, 



(9) 



but as the form (fj[|) is only approximate one should not 
take higher derivatives too seriously. 

Going further down in scale, the form of n(N) depends 
on the assumed interactions that determine the RGE's. 
Typically, n' will increase or decrease monotonically. As 
we have seen, black hole formation will then need n' to 
have a roughly constant value, n' ~ 10~ 2 , and this can 
be achieved with suitable interactions [20( . From Eq. , 
n' ~ 10~ 2 corresponds to to c ~ s ~ ±^/n /2, making 
|c| ~ 10 _1 , in agreement with the expectation ([5]). 



III. RUNNING MASS MODEL 

* 2 See for instance [2ll | and references therein for models with strong 
negative running 

Now we assume that C is generated by the inflaton per- #3 This corrects the re i at ion n' > (n - 1) 2 /4 given in Q. We are 

turbation in a single-field slow-roll model. Of the many ignoring the correction to slow roll invoked in part of that work. 
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IV. OBSERVATIONAL BOUND ON THE 
RUNNING MASS MODEL 

The most recent comparison of the running mass model 
with observation was made in 2004 using WMAP (year 
on) and galaxy survey data available at the time [15| . It 
gave n' < 0.04 or so, easily allowing black hole formation 
in versions of the model where n' does not increase too 
strongly going down in scale. In this section we report 
an update to the earlier bound, using year three WMAP 
data and more recent galaxy survey results. 

In the earlier fit, we took c and s as the parameters 
to be fitted, and only afterward generated contour plots 
of no versus n' . In the present fit, we instead took no 
and n' as the parameters to be fitted. Taking advantage 
of the fact that Eq. ^ practically excludes negative n , 
and that it requires n' to have slow variation, we took n' 
to be constant and imposed n' Q > as a prior. As in the 
previous fit we took the tensor perturbation to be negli- 
gible since that is a prediction of the model. This differs 
substantially from the method adopted in 0] where the 
running of the spectral index was let free to negative val- 
ues and where tensors were included. In that case a neg- 
ative value of the running is obtained, with no running 
excluded at the level of - la (see e.g. @, H d GJ, d ) • 
While this approach is obviously correct when a general 
set of inflationary models is considered, it is important 
to stress that in our case, where we don't consider mod- 
els with n' < 0, the inclusion of those models could bias 
the result towards more restrictive bounds. Moreover, a 
model with n' = gives an acceptable goodness-of-fit to 
the WMAP data and it is therefore statistically legiti- 
mate to assume the prior n' > 0. 

As is now common practice, we base our analysis on 
Markov Chain Monte Carlo methods making use of the 
publicly available cosmomc package f22|. We sample 
the following dimensional set of cosmological parame- 
ters, adopting flat priors on them: the physical baryon 
and CDM densities, ujb — Qbh 2 and uj c = fl c h 2 , the ra- 
tio of the sound horizon to the angular diameter dis- 
tance at decoupling, Sl the scalar spectral index, n, and 
the optical depth to reionization, r. We consider purely 
adiabatic initial conditions. We choose a pivot scale at 
k = 0.002/i- 1 Mpc. 

The MCMC convergence diagnostics are done on 7 
chains applying the Gelman and Rubin "variance of chain 
mean" /"mean of chain variances" R statistic for each 
parameter. Our 1 — D and 2 — D constraints are ob- 
tained after marginalization over the remaining "nui- 
sance" parameters, again using the programs included in 
the cosmomc package. Temperature, cross polarization 
and polarization CMB fluctuations from the WMAP 3 
year data @, [H, [24], [25| are considered and we include a 
top-hat age prior 10 Gyr < t n < 20 Gyr. 

We also consider the small-scale CMB measure- 
ments of the CBI [21, VSA [27|, ACBAR [H and 
BOOMERANG-2k2 [29| experiments. We combine the 
CMB data with the real-space power spectrum of galax- 
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FIG. 1: 68% and 95% c.l. likelihoods in the n - n' plane 
from the WMAP data alone (Top Panel) and WMAP+LSS 
(Bottom Panel). 



ies from the 2dF survey [30j. We restrict the analysis to 
a range of scales over which the fluctuations are assumed 
to be in the linear regime (technically, k < 0.2/1" 1 Mpc) 
and we marginalize over a bias b considered as an addi- 
tional nuisance parameter. 

In Figure [1] we plot the 68% and 98% confidence levels 
in the n-n' plane for two different choices of our datasets: 
WMAP data alone, that should be considered as the 
most conservative result, and "WMAP+LSS" that in- 
cludes the remaining CMB data and galaxy clustering 
data from 2dF. 

As we can see from the Figure, when negative running 
is not considered, the data is still in good agreement with 
a small, but still non-zero running. When the WMAP 
dataset is considered we found a 95% c.l. upper limit 
of n' < 0.039, while the spectral index is bound to be 
n = 0.935±g g|| again at 95%. The best fit (maximum 
likelihood) model has a negligible running n 1 = 0.005 
and n = 0.953. When the remaining cosmological data 
are included, we found a stronger bound on running, with 
n' < 0.026, and n = 0.940lg;^ at 95% c.l.. The best fit 
(maximum likelihood) model has parameters n' = 0.0026 
and n — 0.951. We conclude that the value n' ~ 10 -2 
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creased monotonically as we go down in scale, but black 
hole formation may still be possible with a more compli- 
cated behavior of n'. We are going to exhibit a couple of 
forms of n(k) that would do the job, and still be compat- 
ible with slow- roll inflation. In order to do that, we need 
to consider carefully what the slow-roll approximation 
involves. ^ 4 

The slow roll formalism is reviewed for instance in 
0, [H, HH ■ It starts from the exact Friedmann equa- 
tion #5 



3M 2 H 2 = V{<j>) + \it>\ 
and the exact unperturbed field equation 
4> + 3H<j)+ ^- = 0, 



from which follow the identity 



-0 2 /2- 



(10) 



(11) 



(12) 



FIG. 2: The band corresponds to the spectrum Vq versus In k, 
with constant slope corresponding to n = 0.94. The width of 
the band corresponds to a fractional uncertainty 0.025. The 
range of In k corresponds to the range of cosmological scales, 
explored by observation of the cmb anisotropy and galaxy 
surveys. We see that the band is narrow enough to make the 
variation of Vq significant over the cosmological range. 



required by the running mass model is still viable. 

We have checked that our limits on n' are also con- 
sistent with the WMAP results even in the case when a 
negative running is allowed. The fitted parametes are in 
reasonable agreement with the limits we quoted above, 
even if slightly more stringent due to the inclusion of 
negative running. 

In performing this fit, we chose a pivot point k = 
0.002/i Mpc -1 , corresponding to AN = 1.8. Our fitted 
value n — 0.94 therefore corresponds to no = 0.94 — 1.8n'. 
(Recall that the subscript denotes the scale k = Hq = 
0.00033/i Mpc^ 1 .) In Figure El we plot the shape of 
the spectrum with (i) no running and n = 0.95, (ii), 
n' = 0.01, and n = 0.95 at the pivot point. In the second 
case, n — 1 soon climbs to positive values as we go down 
below the pivot scale. (The finite width of the band is not 
important at this stage, and will be discussed in Section 
VI.) 



In its most basic form, the slow-roll approximation con- 
sists of the two assumptions; 



< 1 



H<t> 



< 1. 



(13) 



These assumptions are usually stated in the equivalent 
forms 



3i70 ~ 

and they imply e<l where 



V{<f>) 
-V'(cf>), 



(14) 
(15) 



(16) 



To calculate the spectrum V((k) one usually considers 
additional parameters; 



9 V" 

V - Ml- 



e = Mp 



a 3 = Mp 



V'V" 

V' 2 V"" 
V 3 ' 



(17) 
(18) 
(19) 



The last two parameters can have either sign despite the 
notation. 



V. SLOW ROLL FORMALISM 

In the future, observation may require negligible run- 
ning on cosmological scales. We have seen that this would 
not permit black hole formation if n' increased or de- 



# 4 We consider only single-field inflation models. The generation 
of black holes has recently been investigated within a double 
inflation model with a strong negative running 31]. 

# 5 As usual Mp = 2 X 10 18 GeV is the reduced Planck mass and V 
is the potential of the inflaton field </>. 
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Using Eq. JT5J) one finds 
din if 



dN 
d(lne) 
d^V 
di] 

dN 



-4e + 2r/ 



£ - -(4 e -,)fW. 



More generally one can define 



am = M P 



2m y / m ~ 1 (d m + 1 y/d^ m + 1 ) 



(20) 
(21) 
(22) 
(23) 

(24) 



which satisfy 
da m 



dN 



[(m - 1)77 - 2me] a m + a m+ i , (25) 



Assuming that the first derivative of Eq. (|15p is also a 
valid approximation, one finds 



(26) 



Comparing Eq. (fl~5|) with the exact equation, we see that 
the fractional error in Eq. (fT5"|) is 0(e, 77), and so we re- 
quire the additional condition 



\t]\ < 1. 



(27) 



Assuming that the curvature perturbation is generated 
from the vacuum fluctuation of Sep, its spectrum in the 
slow roll approximation is given by 35] 



V c {k) 



1 



V 



247T 2 Mp= e 



[l + <D(e,r,)}, 



(28) 



V 3 



12tt 2 M6 



[l + 0(e,77)] (29) 



The right hand side is evaluated at the epoch of hori- 
zon exit. The displayed uncertainty takes account of the 
fractional error in the slow roll approximation that we 
just estimated, and first order (linear) corrections to the 
calculation of the vacuum fluctuation described by the 
Mukhanov-Sasaki equation [3?], [H, [H, H(| (It docs not 
account for nonlinear effects, coming from interactions of 
5cj>. Such effects, which would gen erated non-gaussianity 
of £, are expected to be small [41j . 

Now differentiate Eq. p8|) with respect to Ink, using 
dhik = —dN and ignoring the uncertainty. One finds 

M 



n — 1 = 2i] — 6e 

n' = 2£ 2 + 24e 2 - 16e?7. 



(30) 
(31) 



The error in n — 1, coming from the derivative of the 
error in V^, is 0(£ 2 , e 2 , erf). We will assume that there is 



no cancellation between the two terms of Eq. (|30"|) . and 
make a similar assumption for n' and higher derivatives. 
Also, we will assume that e is negligible compared with 
77, £ 2 and any other relevant a m . Then the fractional 
uncertainty in n — 1 will be small if and only if 



l£ 2 l«M- 



(32) 



Similarly, the fractional error in n' will be small if and 
only if 



«ie 2 



(33) 



In principle one can go on to calculate higher derivatives 
of n, requiring a more extended hierarchy 



\am+i\ < \a m \. 
From Eq. (|2"Tj) . this is equivalent to 



d m lne 



dN* 



L ln« 



dm- 



(34) 



(35) 



We have been exploring the validity of successive 
derivatives with respect to A, of the slow-roll approx- 
imation Eq. (|30|) for n — 1. Barring cancellations, the 
validity of these up to a given order will be equivalent to 
the validity of derivatives of the slow-roll approximation 
for the field equation Eq. (jTTJ) , up to one higher order. To 
see this, start with Eq. (|26| which expresses the validity 
of the first derivative of Eq. |TT|). Put it into Eq. (fPT|) . 
and use ln(l + x) ~ x to find the approximation 



ln(|V'|) = ln(3JT|0| 



n 



(36) 



Assuming that the derivative of the approximation 
Eq. (|2"6"|) is also valid, we can differentiate this with re- 
spect to N: 



d\n\V'\ dln(3ff</>) 



dN 



dN 



d 
dN 



>1 



3 



(37) 



Barring cancellations, the first term on the right hand 
side is of order n — 1 as is easily seen by comparing it 
with the derivative of Eq. (|29p . Therefore, barring can- 
cellations, the validity of this approximation is indeed 
equivalent to the validity of the first derivative of the 
approximation Eq. (|30p . and so on for higher derivatives. 

The equivalence that we saw in the last paragraph 
means that the standard slow-roll approximation for 
n — 1, n' and so on will be valid, if the second, third and 
so on derivatives of the slow-roll approximation Eq. (|26j) 
to the exact field equation Eq. (fTTj) are valid. Reverting 
to our assumption that e is negligible, we conclude that 
the hierarchy Eq. (|32|. Eq. (|33|) etc. will hold (justifying 
the standard formulas for n — 1, n' etc.) to the extent 
that derivatives of the slow-roll approximation Eq. (|2"6"]) 
hold. 

With the hierarchy in place, one can systematically im- 
prove the predictions ([2HD , (SOD and (EU) [H, HI] . The hi- 
erarchy is in general satisfied by the running mass model, 
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and including the leading order correction [38(, the run- 
ning mass prediction ^ becomes [T3| 

n[N) 2 - 1 = (s +l.Q6cs)e« N °- N ) -c (38) 
~ (s + 0.50n') e c{N "- N) - c. (39) 

Such corrections are usually equivalent to a change in 
parameters whose values are not known (in this case, a 
change to s), making them of limited practical impor- 
tance. 

Of course, a given inequality in the hierarchy will 
fail for a few Hubble times if its right hand side passes 
through zero. For instance, Eq. (|3"2"j) will fail if r\ passes 
through zero. Then, if e is negligible, n — 1 will pass 
through zero as well, and while it is doing so the frac- 
tional error in its predicted value will become large. Ac- 
cording to our fit to the data, n(N) — 1 will indeed pass 
through zero on some scale near the bottom end of the 
cosmologically accessible range, if n' has a slowly vary- 
ing value of order 10 -2 . The passage of n — 1 through 
zero need not be a matter of concern, as the absolute er- 
ror remains the same. The running mass prediction (|39[) 
should remain valid even as n passes through zero. 

More generally, it could happen that derivatives of 
Eq. (fT5|) beyond the first are invalid over an extended 
range, so that the hierarchy fails over an extended range. 
To handle such cases one can use the exact (at first order) 
Mukhanov-Sasaki equation or an analytic approximation 

m 



VI. FINITE DIFFERENCE VERSION OF THE 
SPECTRAL INDEX 

Although the hierarchy leads to simple and widely- 
used results, we have seen that its use may sometimes 
be problematic and we will see some more examples of 
that in the following two sections. For that reason, we 
explain in this section how the hierarchy can if necessary 
be avoided. 

The starting point is to realise that the prediction 
Eq. (|2"8")) . with a suitably small error, will accurately de- 
fine the variation of Vq over a finite range, even if the 
mathematical derivative should have large errors (com- 
ing for instance from an oscillation or a break). This is 
illustrated in Figured 

Let us therefore redefine n — 1 so that it specifies a 
finite difference: 



1 



\tiVq2 — In T^i 
Alnfc : 



(40) 



where Alnfc = \x1k2 — hifci and Vqi = V^{ki). To the 
extent that observational bounds on the variation of n(fc) 
are quite weak, this finite difference is really about all 
that observation can determine at present, with Alnfc ~ 
7 or so. 



The error in n — 1 generated by a fractional error x in 
the prediction (j2"8|) will be at most of order 



5{n - 1) 



Alnfc' 



(41) 



Let us assume x ~ 0.025, corresponding to Eq. (|28|) with 
e <tC I77I ~ 0.05 (the observed value of \n — 1|). Then the 
error in the theoretical prediction will be small, if the 
prediction satisfies 



In - II > 



Alnfc 



0.004. 



(42) 



As illustrated in Figure this is very well satisfied if 
h — 1 has the observed value ~ —0.05. 

We can go a bit further, to consider a finite-difference 
version of the running; 

^ f lnTfr - lnP C3 _ InTfe -ln7> ca /(A]nfcjW) 



Alnfc/2 Alnfc/2 
ln^ C2 - 21nP C3 + lnP a ' 
(Alnfc) 2 /4 



(44) 



with 2 In kj, = In k± + In fc2 . The error in n' generated by 
the error x in the prediction will be at most of order^ 6 



Sn 



Wx 



(Alnfc) 2 ' 

The error will be small if the prediction satisfies 
lOx 



W > 



(Alnfc) 2 



5 x 10- 3 , 



(45) 



(46) 



where we again set x = 0.025 as an illustration. Taking 
account of the uncertainties, the prediction for the finite- 
difference version of the running may be valid if n' ~ 
lO" 2 . 



VII. FLOW EQUATIONS 

In the above analysis we worked directly with the po- 
tential. A different approach works initially with the field 
<f>(t), connecting only later with the potential. The start- 
ing point is Eq. (|12[) . providing a parameter e# which 
may be written in various forms; 



H 
1 

2M| 



dilrxH) 

dN 
2 



dN 



, , 1 dH\ 



~jt I >■ (47) 



* 6 The factor 10 accounts roughly for the 1/4 in the denominator 
and the four terms of the numerator. 
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Its derivatives satisfy the exact set of equations 
d(lnejj) 



dN 
dN 



-2(e H + 5i) 
(Si - meH)S m + S m +i, 



where KM k£ 



8 m = d m (j)/dt m /H m 4>. 



(48) 
(49) 

(50) 



Equivalently one can use 4> instead of N as the vari- 
able. Then El 



d(ln(e*) 



d/3„ 



d<p 



-2(e ff -/3i) (51) 
[(to me H ] f3 m + m+1) (52) 



where 142 



A, 



2M F 



dH 



d m+1 H 
d^ 1 ' 



(53) 



These are referred to as flow equations. 

The flow equations (equivalently Eqs. (|4*9"j) and ([50")) ) 
resemble Eqs. ([2D]). ([2"Tjl. and (g5|) but are exact. Slow 
roll with the potential hierarchy ([34[) . up to m = M, is 
obtained if there is a hierarchy* 7 



or equivalently 



|<Wr-l| ^ l^n 



|/3m+l| *C |/3m| 



(54) 



(55) 



up to m = M + 1 . Following [43| one might call this the 
'Hubble hierarchy', as opposed to the 'potential hierar- 
chy' dsn. 

Conversely, if the potential hierarchy is satisfied up to 
m = M, then one can expect the solution 4>(t) to satisfy 
the hierarchy (|55p (equivalently ([Ml)) up to m = M + 1, 
at least for low M. This is because the slow- roll approx- 
imation (|15p is known to be a strong attractor for a wide 
range of initial conditions. As with the potential hierar- 
chy, the Hubble hierarchy will fail briefly if a parameter 
(<5 m or j3 m passes through zero, and it might not be valid 
at all. 



# 7 The stronger hierarchy |5 m +i| 1 /( m + 1 ) < l^ml 1 ^" 1 ' is sometimes 
considered (equivalently |/3 m+ i| 1 /(™.+ i) < |/3 m |V("->). It implies 



the potential hierarchy 



|a m+1 |V(™+i) < 



which is 



satisfied by a wide class of potentials but not by the running 
mass potential. 



VIII. TWO FORMS FOR THE POTENTIAL 

Now we consider forms of the potential, which would 
permit slow roll inflation leading to black hole formation, 
and be consistent with a negligible value of n' . A com- 
mon procedure for generating potentials consistent with 
assigned values of (say) no and n' Q uses the flow equa- 
tions. The equations are numerically integrated with an 
initial hierarchy imposed such as |/3 m+ i|/|/3 m | < 1/5 [II] 
or 1/10 (42J. This procedure is quite complicated, and 
will obviously miss potentials violating the initial hierar- 
chy as discussed at the end of the previous section. # 8 

Our procedure will be to simply specify suitable forms 
for d(lne)/dN). (It resembles the one advocated in [45j.) 
From these the potential can be constructed using 



e{N) = e(7V )exp 



A', 



dine 



dN 



H(N) 

<j>(N ) - 4>(N) 
V(4>) 



N 



N 



H(N ) exp 
/•No 

M P / ^2e(N)dN 

J N 

r2 Tj2/ 



e(N)dN 



3M£ir(0). 



(56) 
(57) 

(58) 
(59) 



The value H(Nq) is determined once the inflation scale 
V(N ) is set, and then e^o) is obtained from Eq. ([28]) 
using the observed value V((N ) = (5 x 10~ 5 ) 2 . 

To keep things simple we focus on small-field inflation, 
which corresponds to e far below 1. To be precise, we 
assume e <C I/Nq for < N < No, corresponding to 
tensor fraction r <C 16/iVo ~ 0.03. Then V can be taken 
to be constant, and black hole formation requires 



Aim 



ln[e(JV )/e(0)] ~ ln[P c (0)/V c (JV )] s 14. (60) 



The predictions for the spectral index and its running 
are then 



n- 1 



= « 



9 V" 



1 rflm 



»' ,2 , r4 V'V" 

2 = ^ ~~V 2 ~ 



2 dN 

Id 2 lne 



2 dN 2 



(61) 
(62) 



We have considered following two forms 
d(lne) „fN\ p f^ N x q 



dN 



B 



Nn 



1 



No 



D : Case I, 



(63) 



* s The flow equations were used in |44| to search for potentials 
consistent with black hole formation but none were found. The 
authors concluded that "... it seems extremely unlikely to us that 
primordial black holes formed as a result of inflationary dynam- 
ics". It was the discrepancy between this result and the earlier 
positive conclusion of 12011 that motivated the present investiga- 
tion. We suppose that it is caused by the use in [44j of the flow 
equations and the hierarchy. 
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n(N ) 


n(N - ATV) 


n'(JVb) 


n'(N - A TV) 


Case I 0.9500 


0.9529 





0.0017 


Case II 0.9500 


0.9511 


1.4942 x 10" 15 


5.9700 x 10" 16 



TABLE I: n and n' at TV 
and Case II. 

and 



TV and TV = TV - ATV in Case I 



ef(lne) 
dN 



B exp 



TV 



N — ATV - A 



D : Case II. 



(64) 



In Fig.[3]we plot the schematic pictures in case I (top) and 
case II (bottom) , respectively. In figure [5] we show the 
derivatives of d In e/dN with respect to TV. The hierarchy 
(|35p is in general respected except where the denomina- 
tors pass through zero. 

We impose the observational constraints on Vq, n and 
n' . To be on the safe side we also impose a rough finite- 
difference version of the constraint on n' in the following 
way. The WMAP data spans a range roughly £ ~ 0(1) - 
0(1O 3 ), corresponding to ATV <~ 7 (= ATV, and n derived 
from that data has an error of about 0.1. Therefore, 
we require that n should change by less than 0.02 in 
the range TVo to TVo — ATV. We have checked that the 
condition (|60jl needed for the PBH production is satisfied 
with all of the observational constrains for p = 1, q = 3, 
B ~ 5.5, D ~ 0.05, TV = 60 and ATV = 5 for Case 
I, and q = 10, A = 5, B ~ 0.5, D ~ 0.05, TV = 50 
and ATV = 10 for Case II. For the detailed numerical 
values of n and n', see Table |TJ It is clear that with a 
parameterisation like the one in Figure [31 we can make 
n practically constant over the range ATV ~ 10 of scales 
probed by large-scale observations. 

For the parameterisation I, the potential has a scal- 
ing V cx (<f> - <j) Q ) 2 . We plot it in Figure |U Note that 
Vq 1 ^ 4 /M p should be less than ~ 10 -3 for the slow-roll in- 
flation e <C 1/TVo and <f> < M p . This shape is similar to 
those in some classes of hilltop inflation models [46[ . 



IX. BLACK HOLE FORMATION IN A 
CURVATON-TYPE PARADIGM 

The inflaton contribution Ci(^) is time-independent, 
and is the only one present at horizon exit. Subsequently 
though, the contribution C c (£, k) of some other (curvaton- 
type) field could grow and become dominant. (See [8( for 
a discussion of the possibilities with extensive references.) 
Eventually C,{t) will level out to the observed value; 



C(fc) = ( t (k) + Uk), 



(65) 



where the last term is the eventual time-independent 
value of the curvaton-type contribution. In an obvious 
notation, the observed spectrum is now 



(a) Case I 




FIG. 3: Schematic pictures of functional form of y(N) in Case 
I (top) and Case II (bottom), respectively. For reference we 
also plot the constant case, n = 0.3 (see the text). 



N=100 



0.00014 
0.00012 
0.0001 
0.00008 
0.00006 
0.00004 
0.00002 





FIG. 4: Form of the potential V as a function of the field 
d>. The horizontal axis is the normalized value of d>, d> = 



\ 4>-4>a 1 / 
[ M v J / 



with cf> = 0(0) and V = V(<f> ). The 



vertical axis means V 



ft- 1 ]/ 



-2 /10" 



Vt = Vi + v c , 



(66) 
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(a) Case I 




(b) Case II 



o.oi 




FIG. 5: Derivatives of d(\ne)/dN with respect to TV for Case 
I (top) and Case II (bottom), respectively. Here we plot the 
higher derivatives (In e) w = d {l) '(he) / 'dN w for I = 2, 3, and 
4. 




In H„ 



Ink 



In 10H„ 



FIG. 6: Two scenarios for generating black holes using a 
curvaton-type paradigm. 



and the spectral index is 

n-l = fi(rn - 1) + f c (n c - 1), 



(67) 



where fi = VijV^ and f c = Vc/V^. We will consider 
two different possibilities for the ratio Vi(k) /V c {k), illus- 
trated in Figure [Sj 



A. Black holes from the inflaton perturbation 

We first assume that f c -C 1 at the end of inflation so 
that black holes are generated by the inflaton perturba- 
tion, but that /jCl while cosmological scales leave the 
horizon. To agree with observation, we will demand at 

N = N a 



and 



l~7i- 1 ~ 0.05, 



film - 1| < i x lo- 



To form black holes, Eq. (|60|) becomes 



14 = ln[P C (0)/P c (iVo)] 

= ]n[PdO)/n(No)}+Hfi(0)] 



(68) 



(69) 



(70) 
(71) 
(72) 



where we set = to get the last line. These require- 
ments are satisfied with, for example, /i(0) ~ 10 _1 and 
rii ~ 1.4, and the observational bound on the running 
imposes no further constraint. 

The required spectral index rij corresponds to t] = 0.2 
which is more or less compatible with the slow-roll re- 
quirement 1 77 1 <c 1. In the context of supergravity such a 
value is more natural than the small value r\ = —0.025 re- 
quired to fit observation. This looks promising for black 
hole formation, but we have to remember that in a cur- 
vaton type model the prediction for n becomes Q 



Vera 



2r} aa - 2e 
1 d 2 V 



(73) 
(74) 



If e is negligible as we are assuming, this requires r\ aa — 
—0.025 which may difficult to achieve since a will tend 
to roll away from any maximum of its potential [8j . In a 
curvaton type model it may therefore be more attractive 
[H| to take ri aa negligible and e = 0.025, but we have not 
explored that option. 



B. Black holes from the curvaton-type contribution 

Now we suppose that the roles of the inflaton and 
the curvaton are reversed, so that the inflaton gener- 
ates the observation curvature perturbation while the 
curvaton perturbation generates black holes. In this case 
black hole generation occurs only when the curvaton-type 
mechanism operates which will usually be long after in- 
flation is over. 

In this scenario, we have to modify Eqs. (f6"9")) - (|T2"|) by 
interchanging i and c, and replacing the epoch N = by 
an earlier epoch N cnrv . This is the epoch at which the 
scale leaves the horizon, that is entering the horizon when 
the curvaton mechanism operates. To achieve black hole 
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formation we will therefore generally need n c > 1.4, but 
that need not be a problem. Indeed, within the context 
of supergravity a value n c significantly bigger than 1 is 
expected Q, just as it is for n,. This, our third scenario 
for generating black holes, therefore seems at least as 
good as the other two. 

A particularly interesting possibility in this case, is 
that the curvature perturbation generating the black 
holes could easily be highly non-gaussian, to be pre- 
cise the square of a gaussian quantity Q. This would 
not make much of a change [47] in the magnitude of Vq 
needed to generate black holes (upon which our estimates 
are based) but it could alter the predicted shape of the 
black hole mass function. 



X. CONCLUSION 

The possibility of primordial black hole formation at 
the end of inflation has a long history, which was long 
overdue for an update. The update is needed partly be- 
cause observation now requires on cosmological scales a 
tilt far below 0.3 (and with negative sign) and not too 
much running. It is also needed because the original 
paradigm, that the inflaton perturbation is entirely re- 
sponsible for the curvature perturbation, is now only one 
possibility. 

According to the standard paradigm, the curvature 



perturbation is generated during slow roll inflation from 
the vacuum fluctuation of the inflaton field. Within 
this paradigm, the running-mass model provides a well- 
motivated way of achieving black hole formation. To 
form black holes, the model probably requires strong run- 
ning on cosmological scales, n' ~ 0.01. 

Such running is allowed by the data. If a value 
n' ~ 10 ~ 2 on cosmological scales is ruled out in the fu- 
ture, n' will still have to increase to > 10~ 2 in order 
to form black holes. We saw that this may be achieved 
within the standard paradigm by a suitable potential. 
Alternatively, it might be achieved by a switch from the 
standard paradigm to a curvaton-type paradigm, or by a 
switch from the curvaton-type paradigm to the standard 
paradigm. 

For observers, we would like to re-iterate and earlier 
conclusion [3J|, that an upper bound n! < 10 -3 , or de- 
tection, would have important implications and is a very 
worthwhile goal. 
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